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ABSTRACT 

It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be 
discovered. Water content helps determine surface land fraction, which influences planetary weathering 
behavior. This is important because the silicate weathering feedback determines the width of the 
habitable zone in space and time. Here a low-order model of weathering and climate, useful for 
gaining qualitative understanding, is developed to examine climate evolution for planets of various 
land-ocean fractions. It is pointed out that, if scafloor weathering does not depend directly on surface 
temperature, there can be no weathering-climate feedback on a watcrworld. This would dramatically 
narrow the habitable zone of a waterworld. Results from our model indicate that weathering behavior 
does not depend strongly on land fraction for partially ocean-covered planets. This is powerful because 
it suggests that previous habitable zone theory is robust to changes in land fraction, as long as there 
is some land. Finally, a mechanism is proposed for a waterworld to prevent complete water loss 
during a moist greenhouse through rapid weathering of exposed continents. This process is named 
a "waterworld self-arrest," and it implies that waterworlds can go through a moist greenhouse stage 
and end up as planets like Earth with partial ocean coverage. This work stresses the importance of 
surface and geologic effects, in addition to the usual incident stellar flux, for habitability. 
Subject headings: astrobiology - planets and satellites: general 



1. INTRODUCTION 

The habitable zone is traditionally defined as the re- 
gion around a star where liquid water can exist at the 
surface of a planet (jKasting et al.l fl993f ). Since climate 
systems include both positive and negative feedbacks, a 
planet's surface temperature is non-trivially related to 
incident stellar flux. The inner edge of the habitable 
zone is defined by the "moist greenhouse," which oc- 
curs if a planet becomes hot enough (surface tempera- 
ture of «340 K) that large amounts of water can be lost 
through disassociation by photolysis in the stratosphere 
and hydrogen escape to space. The outer edge occurs 
when CO2 reaches a high enough pressure that it can 
no longer provide warming, either because of increased 
Rayleigh scattering or because it condenses at the sur- 
face, which results in permanent global glaciation. These 
limits do not necessarily represent hard limits on life 
of all types. For example, life survived glaciations that 
may have been global ( "Snowball Earth" ) on Earth 600- 
700 million years ago (|Kirschvinld IT992: Ho ffman et al.l 
Il998f ). Furthermore, if greenho use gases other than CO2 
are c onsidered, e.g . , hyd rogen (jPierrehumbert fe Gaidosl 
120111: IWordswortbl 120121 ), the outer limit of the habit- 
able zone could be pushed beyond the CO2 condensation 
limit. Alternative limits on habitability have also been 
proposed. For example, a reduction in the partial pres- 
sure of atmospheric CO2 below 10 -5 bar may prevent C4 
photosynthesis, which could curb the complex biosphere 
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(jCaldeira fe Kastind 119921: iBloh et alll2005l) . although it 
would certainly not limit most types of life. 

Studies of planetary habitability have become in- 
creasingly pertinent as capabilities for the discovery 
and characterization of terrestrial or possibly terres- 
trial exoplanets in or near the habita ble zone have im- 
proved. For exami al e, planets GJ581 d (lUdrv et alll2007l: 
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(| Anglada-Escudi et all 1201 2h have recently been discov- 
ered and determined to be potentially in the habitable 
zone. It is likely that the number of good candidates for 
habitable terrestrial exoplanets will increase in the near 
future as the Kepler mission, ground-based transit sur- 
veys, radial velocity surveys, and gravitational microlens- 
ing studies continue to detect new planetary systems. 

The "Faint Young Sun Problem" discussed in the 
Earth Science literature is analogous to the habitable 
zone concept. A planet orbiting a main-sequence star re- 
ceives ever increasing incident flux ("insolation") as the 
star ages. This should tend to hurry a planet through the 
habitable zone unless the planetary albedo or thermal op- 
tical thickness (greenh ouse effect) adjus t as th e system 
ages. As pointed out bv lSagan fc Mullenl (|1972D . geologi- 
cal evidence for liquid water early in Earth's history sug- 
gests that just such an adjustment must have occurred 
on E arth. Although there is some debate about the de- 
tails (|Kastingll2010( ). it is fair ly widely accepted t hat the 
silicate- weathering feedback ([Walker et al.lfl98ll ) played 
an important role in m aintaining cle ment conditions 
through Earth's history l)Feulnerl 120121 ). Through this 
negative feedback the temperature-dependent weather- 
ing of silicate rocks on continents, which represents the 
main removal process for CO2 from the atmosphere, is re- 
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duced when the temperature decreases, creating a buffer- 
ing on changes in temperature since CO2 is a strong in- 
frared absorber. The silicate weathering feedback also 
greatly e xpands the habitable zone annulus around a star 
m space (|Kasting et al.lll993D . since moving a planet fur- 
ther from the star decreases the insolation it receives just 
as moving backward in time does. If insolation bound- 
aries are used to demarcate the habitable zone, the con- 
cept can be used to consider habitability as a function of 
position relative to the star at a particular time, or as a 
function of time for a planet at a constant distance from 
a star. 

In addition to continental silicate weathering, weath- 
ering can also occur in hydrothermal systems in the 
basaltic oceanic crust at the seafloor. Scafloor weather- 
ing is poorly constrained on Earth; however, it is thought 
to be weaker than continental weathering and to de- 
pend mainly on ocean chemistry, pH, and circulation 
of seawater throug h basaltic crus t , rather than directly 
on surface climate (|Caldeiralll99l iSleep fe ZaimTell2001l; 
iLe Hir et all I2008D . If this is true, C0 2 would be less 
efficiently removed from the atmosphere of a planet with 
a lower land fraction, leading to higher CO2 levels and 
a warmer climate. Furthermore, a planet with a lower 
land fraction could have a weaker buffering to changes 
in insolation than a planet with a higher land fraction, 
which would cause it to have a smaller habitable zone. 

While the weathering behavior of a planet could de- 
pend on surface land fraction, the water complement of 
planets in the habitable zone should vary substantially. 
The reason for this is that the habitable zone is in gen- 
eral located closer to the star than the snow line, the 
location within a protoplanetary disk outside of which 
water ice would be present and available to be incorpo- 
rated into solids. In the solar system, for example, the 
current habitable zone rang es from approximately 0.8 to 
1.7 AU (IKasting et al.lH99l if one neglects the uncertain 
effects of CO2 clouds. If CO2 clouds are assumed to pro - 
duce a strong warming (|Forget fc Pierrehumbertlli997l ). 
the outer limit of the habitab l e zone could be extended 
to »2.4 AU (IMischna et alJ 120001 ISelsis et all [2001 
IKaltenegger fc Sasselovl 120 111) , although recent three- 
dimensional simulations including atmospheric CO2 con- 
densation suggest that the warming effect may in fact 
be fairly modest dWord sworth et al.l 120111 : iForget et al.l 
IMllWordsworth et al.1120121) . The snow line is though t 
to have been around ~2.5 AU (jMorbidelli et al.ll2000T ); 
beyond this distance, evidence for water in the form of 
ice or hydratcd minerals is seen in asteroids and plan- 
ets. The general picture for water delivery to Earth has 
been that the orbits of bodies from beyond the snow line 
were excited to high eccentricities through gravitational 
scattering by massive planetary embryos and a young 
Jupiter. These high eccentricities would have put water- 
bearing planetesimals and embryos on Earth-crossing or- 
bits. A fraction of such bodies would have been ac- 
creted by Ear th, stochastically delivering volatiles to the 
young planet (IMorbidelli et al.ll2lMlO'Brien et"ail[200l 
IRavmond et al.H2009D ~ 

Efforts have been made to extend this analysis to plan- 
etary systems around other stars. The same dynamical 
effects are expected to occur, with the amount of water 
delivered to a potentially habitable planet being strongly 
dependent on the presence or orbits of giant plan- 



ets dO'Brien et al.ll200dlRavmondll2006t IRavmond et al.1 
120091). As low mass stars are expected to have lower 
mass disks and therefore fewer massive bodies to pro- 
duce gravitational scattering, low mass stars may be 
more likely to have vo latile-poor habitable-zone planets 
([Raymond et al.ll2007f ). Stars of mass ~1 M Q or higher, 
on the other hand, could easily have habitable-zone plan- 
ets of similar or greater fra ctional mass of water than 
Earth (|Ravmond et al.ll2007t ). 

As described above, current thinking mostly focuses 
on hydrated asteroids as the main source of water for 
a habitable planet; however, there are other dynamical 
mechanisms that could allow habitable planets to accrete 
significant amounts of water and volatiles. Comctary- 
bodies would deliver a larger amount of water for a given 
impactor mass than asteroids. Should such bodies get 
scattered over the course of planet formation, this may 
allow larger amounts of water t o be accr e ted th an oth- 
erwise predicted. Furthermore, iKuchnerl (|2003t ) argued 
that planets which formed outside the snow line could 
migrate inwards due to gravitational torques from a pro- 
toplanetary disk or scattering with other planets. Such 
bodies would naturally accrete large fractions of water 
ice during formation beyond the snow line, which would 
then melt and sublimate in a warmer orbit, providing a 
high volatile content to a planet close to its star. 

Habitable zone planets could therefore have a wide va- 
riety of water mass fractions, which would lead to varying 
land fractions. Furthermore, even for a constant mass 
fraction of water, scaling relations dictate that land frac- 
tion depends on planetary size. At the same time, the 
surface land fraction could exert strong control on the 
planetary carbon cycle, which strongly influences plane- 
tary habitability. This warrants a general consideration 
of weathering behavior on planets of varying land frac- 
tion. 

The main objective of this paper is to investigate the 
effect of land fraction on the carbon cycle and weathering 
behavior of a terrestrial planet in the habitable zone. We 
will outline and use a simple analytical model for weath- 
ering and global climate that necessarily makes grave 
approximations to the real physical processes. For ex- 
ample, we will use existing parameterizations of scafloor 
weathering, while acknowledging that observational and 
experimental constraints on such parameterizations are 
minimal. We will only use the model, however, to make 
statements that do not depend strongly on uncertain as- 
pects of the parameterizations. This model should be 
used to understand intuitively the qualitative behavior of 
the system rather than to make quantitative estimates. 
A major strength of the model is that it is easy to derive 
and understand, yet should capture the most significant 
physical processes. This type of modeling is appropriate 
in the study of exoplanets, for which limited data that 
would be relevant for a geochemical model exist. 

We will consider an Earth-like planet with silicate 
rocks, a large reservoir of carbon in carbonate rocks, and 
at least some surface ocean. We will use equilibrium re- 
lations for weathering, which is reasonable for the slow 
changes in insolation that a main-sequence star experi- 
ences. The climate model we use is a linearization of 
a zero dimensional model. Although this is a severe ap- 
proximation, it allows the analytical progress that we feel 
is useful for obtaining insight into the problem. We will 
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consider planetary surface land fractions ranging from 
partial ocean coverage to complete waterworlds. In this 
context a planet would be a waterworld if the highest 
land were covered by even 1 m of water, although a planet 
with more water than this would qualify as a waterworld 
as well. However we will assume the geophysical context 
is a planet with a substantial rock mantle and only up to 
rough ly »10 times Earth's water ma ss fraction (0.02%- 
0.1%; [Hirschmann fc Kohlstc dt 2012) rather than poten- 
tial w aterworlds that might be 0(10%) or more water by 
mass (|Fu et al.ll2010l ) and would therefore have vastly dif- 
ferent volatile cycles. Although partially ocean-covered 
planets have rec ently been referred to as "aquaplanets" 
(|Abe et aTll2011f ). we will not adopt this terminology be- 
cause "aquaplanet" has a long history of being used to 
denote a completely ocean-covered planet in the climate 
and atmospheric dynamics communities. 

We will find that the weathering behavior is fairly in- 
sensitive to land fraction when there is partial ocean 
coverage. For example, we will find that weathering 
feedbacks function similarly, yielding a habitable zone 
of similar width, if a planet has a land fraction of 0.3 
(like modern Earth) or 0.01 (equivalent to the combined 
size of Greenland and Mexico). In contrast, we will find 
that the weathering behavior of a waterworld is drasti- 
cally different from a planet with partial ocean coverage. 
If scafloor weathering depends mainly on ocean acidity, 
rather than planetary surface temperature, no weath- 
ering feedback operates on a waterworld and it should 
have a narrow habitable zone and progress through it 
quickly as its star ages. Finally, we will argue that it 
is possible for a waterworld to stop a moist greenhouse 
in progress when continent is exposed by drawing down 
the CO2 through massive weathering, which would leave 
the planet in a clement state with partial ocean cover- 
age. We will refer to this possibility as a "waterworld 
self-arrest." 

This paper comp lements in two ways the recent work 
of lAbe et alj (|201lD . who found that a nearly dry planet 
should have a wider habitable zone than a p lanet with 
some water. First, the calculations made by lAbe et all 
(|2011l ) focused on climate modeling rather than weath- 
ering, although they included a qualitative description of 
factors that would influence weathering on a dry planet. 
Second, we consider a variety of land fractions up to the 
limiting case of a waterworld. 

The outline of this paper is as follows. We describe 
our model in Section [2j use it in Section [3j and perform 
a sensitivity analysis in Section |U We discuss the possi- 
bility of a waterworld self-arrest in Section [5] We outline 
observational implications and prospects for confirma- 
tion or falsification of our work in Section [6l We discuss 
our results further in Scction[7J including considering the 
limitations of our various assumptions, and conclude in 
Section [5] 

2. MODEL DESCRIPTION 

Here we will briefly review silicate weathering and the 
carbon cycle before developing our mo del. The read er in- 
terested in more detail shou l d cons ult iBernerl (|2004l ) and 
Chapter 8 of lPierrehumbert) (|2010f) . The carbon cycle on 
a terrestrial planet can be described as a long-term bal- 
ance between volcanic outgassing of CO2 and the burial 
of carbonates, mediated by chemical weathering of sili- 



cates. Silicate minerals are weathered through reactions 
such as 

CaSi0 3 (s) + C0 2 (g) ^ CaC0 3 (s) + Si0 2 (s), (1) 

where we have used CaSiOa as an example silicate min- 
eral, but MgSiC>3, FeSiOa, and much more complex min- 
erals can also participate in similar weathering reactions. 
The reaction product CaC03 is an example carbonate 
and silicon dioxide (Si0 2 ) is often called silica. These 
reactions occur in aqueous solution and lead to a net 
decrease in atmospheric CO2 if the resulting CaC03 is 
eventually buried in ocean sediment, ultimately to be 
subducted into the mantle. CaC0 3 is generally produced 
in the ocean by biological precipitation, but if there were 
no biology weathering fluxes into the ocean would drive 
up carbonate saturation until abiotic CaC03 precipita- 
tion were possible. Carbon in the mantle can be released 
from carbonate minerals if it reaches high enough tem- 
peratures, leading to volcanic outgassing and closing the 
carbon cycle. The mantle reservoir of carbon is large 
enough that C0 2 outgassing does not depend on the 
amount of carbon subducted into the mantle. 

Weathering reactions can occur either on continents 
or at the seafloor and the weathering rate is the to- 
tal amount of C0 2 per year that is converted into car- 
bonate by reactions like Equation ([1]) and subducted. 
The weathering rate is often measured in units of kilo- 
grams carbon per year (kg C yr _1 ). Below we will work 
with the dimcnsionlcss weathering rate, which is nor- 
malized such that it equals unity at the modern weath- 
ering rate. Because weathering reactions are aqueous 
and because rain increases erosion, allowing more rock 
to react, the continental weathering rate should increase 
with increasing precipitation. Experiments suggest that 
weathering reactions should increase with increasing C0 2 
concentration, but that this dependence may be signif- 
icant l y reduced when land plants are present (|Bernerl 
l2004t iPierrehumbertl I2010T ) . Finally, weathering rates 
are generally assumed to have an exponential depen- 
dence on temperature . This is based on temperature fits 
([Marshall etaD 11988ft to cation. e.g .. Ca 2+ and Mg 2+ , 
data in river runoff (|Mevbecklll979l f. and appears to be 
ultimately du e to an Arrhenius law for rock dissolution 
(jBemerlll994h . 

These dependencies of continental weathering can be 
combined into the following approximation for Wi, the 
continental silicate weathering rate, 




(2) 



where P is the precipitation rate, cf> is the partial pres- 
sure of C0 2 , T is the surface temperature (see Table [1] 
for a list of important model variables), and T u = 10° C, 
a = 0.65, and b = 0.5 are constants. The subscript 
represents the current value and a tilde represents a di- 
mensionful quantity. The precipitation dependence of 
continental weathering, a, is determined experimentally 
by measuring the Ca 2+ and Mg 2+ solute conce ntration 
in riv ers an d regressing against a nnual runoff. IBernerl 
(|1994D and IPierrehumbertl (|2010j) use a = 0.65 based 
on a study in Kenya that found a = 0.6 (|Dunnd[l978l) 
and a study in the United States that found a = 0.7 
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TABLE 1 

List of the Most Important Model Variables 
Variables Definition 

T Global mean surface temperature 

<j> Partial pressure CO2 in the atmosphere 

S Stellar flux at planetary distance 

(solar constant for Earth) 
Planetary land fraction 
Climate sensitivity to external forcing 
Weathering rate 
r Outgassing rate 

(|PeterdH98l . We will use a = 0.65 as our standard 
value, but we find in Section 3] that our main results 
are robust for a between and 2. The C02-dependence 
of continental weathering, b, could be as low as if 
land plants eliminate th e effect of CO2 on weather- 
ing (jPierrehumbertl l2010f) . an d the maximum val ue it 
should take is 1 (lBernerl ll994). Both (H erne rl [T99l and 
iPierrehumbertl (|2010D adopt the intermediate value of 
0.5, as we will here. Our main results are only sensi- 
tive to b if it becomes very small (Section |4} . Various 
estimates of the e-folding temper ature, T,,,, exist in th e 
literature, in cluding T u = 10.0 K (lPierrehumbertll2010h, 
T u = 11.1 K (lBerner| [l99l. T„ = 21.7 K ([Marshall et al.l 
[l988lh and T u = 17.7 K (|Walker et al.lll981lh We adopt 
T u = 10.0 K and show that our main conclusions are 
valid for T u = 5 - 25 K in Section H 

Equation ([2]) can be nondimensionalized to take the 
following form 

Wi = P a /e T , (3) 



where Wi 



Wi 

w l0 ■ 



p_ 

Pa 



and T = ^=^-. 



We parameterize the precipitation as a function of tem- 
perature as 



£ = l + 5p(f -T ), 



(4) 



where a p represents the fractional increase in precipi- 
tation per degree of warming. Climate change simula- 
tions in a variety of global climate models suggest that 
dtp = 0.02 — 0.03, corresponding to a 2%-3% increase in 
global mean precipitation per Kelvin increase in global 
me an temperature, is an appropriate value when T « 
To ([Schneider et al.l \2 01 Of ). Simulations in an idealized 
global climate model, however, suggest that precipitation 
asymptotically approaches an energetically determined 
constant limit at T w 300 — 31 K as the atmospheric 
optic al thickness is increased (|Q 'Gorman fc Schneider! 
•_»()()•>;. The energetic constraint on precipitation would 
be expected to change, however, if the insolation changes, 
which is the situation we will be considering below. For 
simplicity, we will take a p = 0.025 in what follows. 
Because of the exponential dependence of continental 
weathering on T in Equation ([3]), our results are only 
minimally affected by a p , such that our results are not 
substantially altered even if we eliminate the tempera- 
ture dependence of precipitation by setting a p = (Sec- 
tion [4]). Defining ov, = 6t p T u = 0.25 we can nondimcn- 
sionalize Equation (g]) to find 

P = 1 + a p T. (5) 

In addition to continental silicate weathering, reactions 
similar to Equation ([T]), which lead to the net consump- 



tion of CO2, can occur at the scafloor. We will refer 
to this "low-temperature" alteration of basaltic ocean 
crust as "scafloor weathering." We do not consider 
"high-tcmpcrature" basaltic alteration here because it 
does not lead to net CO2 consumption. Scafloor weath- 
ering should depend on many factors, including ocean 
pH, the chemical composition of the ocean, the temper- 
ature at which the reactions occur, circulation of sea- 
water through basaltic crust, and the seafloor spreading 
rate. We will consider a constant seafloor spreading rate 
and assume a chemical composition similar to the mod- 
ern oc ean. Based on laboratory measurements, ICaldeiral 
(|1995l) parameterized the seafloor weathering rate (W ) 
as 

W [H+] m + v 



W o0 [H+]™ 



(O 



where [H + ] is the hydrogen ion concentration and m is a 
constant, v has a value of 10~ 3 and units of M m , where 
M is molar concentration, in order to satisfy dimensional 
constraints in Equation ([5]). This parameterization as- 
sumes that there is enough weatherable material at the 
scafloor, and sufficient circulation of seawater through it, 
to accommodate arbitrary increases in ocean acidity. If 
the weatherable material produced at ocean spreading 
regions is completely used up before subduction occurs, 
W will be limited in a way not included in Equation ((BJ) . 
We will return to this assumption in Section [7] 

Ocean chemistry leads to a functional relationship be- 
tween <f>, the atmospheric CO2 concentration, and [H + ] 
of the ocean. We will assume that the ocean chem- 
istry can be described as a carbonate ion system cou- 
pled to an equilibrium relation for CaC03. Although 
this assumption neglects other sources of ocean alkalin- 
ity, e.g., Mg 2+ , Na + , and K + , it has been useful for mod- 
eling the carbon cycle in deep time Earth problems (e.g., 
iHiggins fc Schr ag 2003). Furthermore, we implicitly as- 
sume that Ca 2+ can be supplied to the ocean by sub- 
marine rocks in a waterworld or a planet with very low 
land fraction. This system is described by the equilib- 
rium relations for the dissolution of CO2 and CaC03, 
the disassociation of H2CO3, HC07, and H2O, and al- 
kalinity (charge) balance (see chapter 4 of iDreverl 119881 
for details). Using the equilibrium constants given by 
IDreverl (fl988h . we find that for an atmospheric CO2 con- 
centration between 10 -4 and 10 2 bar, the dominant al- 
kalinity balance is 2[Ca 2+ ] = [HCO^] in this system. 

Since in the carbonate ion system [Ca 2+ ] oc 



and [HCO~ 3 



pTFp 



[HCO3-] 

this dominant balance leads to 



[H + ] = c(T)(j)3 with c(T) given by a combination of equi- 
librium constants. c(T) varies by only O(10%) between 
0°C and 90°C, and we therefore neglect this temperature 
dependence below. This leads to the nondimensional ap- 
proximation [H + ] = 03 , which allows Equation (j6|) to be 
written in nondimensional form 



W = 



where W Q 



w 

Wor, 



and v = 



the ICaldeiral (|1995ft values (m 



(7) 



Using 



and R = 10" 3 Mi) 
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and assuming a reference ocean pH (pH=-log[H + ]) of 8.2 
yields n = ± and j/=12.6. Note that ICaldeiral (|1995[ ) 
discussed the possibility that an exponent to [H + ] other 
than i may be appropriate and we will discuss other 
scalings in Sections H] and [7] We note that we have ne- 
glected any dependence of seafloor weath ering on surface 
temp erature, which is equivalent (e.g., ISleep &: Zahnld 
to assuming the temperature at which seafloor 
weathering occurs (ps20-40°C) is set by heat flux from 
Earth's interior, rather than by the temperature of sea- 
water. Since this assumption eliminates the potential for 
a weathering-climate feedback due to seafloor weather- 
ing (Section [3]), it represents a conservative assumption 
when considering the effect of changing the land frac- 
tion on weathering. We discuss the effect of allowing a 
seafloor weathering dependence on surface temperature 
in Section [3 

Combining Eqs. ([3]), and ([7]), and assuming that 
continental weathering scales linearly with land fraction 
(7) and seafloor weathering scales linearly with ocean 
fraction (I-7), the total weathering rate can be written 



W = Po-L(l + OpT) Ve J + (1 - A)) 
7o 



1-7 (<l> n + 



l-7o 



(8) 

where 70 is the modern Earth land fraction and /3q = 
w W +w - is the fraction of total weathering accomplished 
on continents for modern Earth conditions (7 = 70, 
T = 0, and 0=1). (3o will not generally be the ac- 
tual fraction of weathering accomplished on continents 
for other climatic conditions. The magnitude of the 
seafloor weathering rate on modern Earth is constrained 
by drilling cores in the ocean floor, measuring the amount 
of basaltic rock that has been weathered, and scaling to 
the entire ocean. Because relatively few core holes have 
been drilled and analyzed, and because weathering rates 
differ significantly among them, the sea floor weatherin g 
rate on Earth is only roughly kno wn (|Caldeira I1995D . 
ICaldeiral (Il99l assumes /3 = 0.67. iLe Hir eTaLl ([2008h 
assume (3q = 0.78, and we will take (3o = 0.75 as our stan- 
dard value, although we show that our results are robust 
as long as f3o ~ 0.5 — 1 in Section|U Our assumption that 
continental and seafloor weathering scale linearly with 
land and ocean fraction, respectively, cannot be rigor- 
ously justified and should be viewed as a first step. It is 
possible, for example, that most seafloor weathering oc- 
curs near spreading centers so that total seafloor weath- 
ering is nearly independent of ocean fraction. If true this 
would not affect our main results significantly since they 
do not depend on the details of the seafloor weathering 
parameterization (Section [4j , but would affect an eval- 
uation of the effect of seafloor weathering on the outer 
edge of the habitable zone, where strong seafloor weath- 
ering could keep CO2 below the limit imposed by CO2 
condensation or Rayleigh scattering (Section [7]) . 
Global mean energy balance for a planet can be written 



S 



(l-a) = n(T, 



(9) 



where S is the dimensionful insolation (50=1365 W m 2 
for modern Earth), a is the planetary Bond albedo, or 
reflectivity («0.3 for modern Earth), 4 is a geometric 



factor representing the ratio of the area of a planet inter- 
cepting stellar radiation to the total planetary area, and 
fi(T, </>) is the infrared radiation emitted to space (out- 
going longwave radiation). Below we will take a — a a to 
be a constant for simplicity. We discuss this assumption 
further in Section [7] 

After nondimcnsionalizing (f2 = j}-, S = ■§-), the emit- 
ted thermal radiation, f2, can be Taylor expanded around 
current Earth conditions (S = 1, T = 0, and <f> = 1) to 
find 

fl(T,cf>)^hl-a)+u} 1 T-uj 2 log^ (10) 



where ui\ 
an 



on 1 

ay|r=o,0=i,s=i 



and uj 2 = 

3i og (0) ir=o,0=i,s=i- Note that for a planet with 
no atmosphere, il oc T 4 , as one would expect for a 
blackbody; the current parameterization accounts for 
the effects of the atmosphere, but is only a linearization 
of more complete radiative transfer. Because the CO2 
absorption bands are saturated, t he outgoing longwav e 
radiation is roughly linear in log 4> (|Pierrehumbertll2010D . 
so wc expand in \og4>. We calculate oj\ and uj 2 using 
published polynomial fits to f2 calculate d using a single 
column radiative-convective model (|Pierrehumbertl 
120101 ). We find wi=0.015 using lji = and calcu- 

lating u>i as the average slope of ft between T = 275 C 
and T — 325 K with CO 2 =1000 ppm and a relative 
humidity of 50%. We find w 2 =0.007 using uj 2 = 
and calculating lJ 2 as the negative of the average slope 
of f2 with respect to \og(j> between CO2=280 ppm and 
C02=2.8 x 10 5 ppm with a surface temperature of 
T=280 K and a relative humidity of 50%. Wc show 
that our main conclusions are unaltered if U3\ and ui 2 are 
changed by a factor or two in either direction in Section 
01 The linearization of f2(T, 0) leads to the following 
climate model (nondimensional form of Equation (j9))) 



1 



-AS - 1)(1 - a) = ujiT - uj 2 log< 



(11) 



Taking T to be the ratio of the CO2 outgassing rate 
to that on modern Earth, CO2 equilibrium in the at- 
mosphere requires V = W , where W is given by Equa- 
tion ((SJ) . Such an equilibrium is set up on a million year 
timescale, so that kinetic (non-equilibrium) effects can 
be neglected when considering variations in the weather- 
ing rate in response to changes in stellar luminosity. We 
make order-of-magnitude estimates on a kinetic process 
in Section [5] and we discuss the possibility that weather- 
ing cannot increase enough to establish r = W in Section 
[7J T = W, along with Equation (fTTj). represents a cou- 
pled climate-weathering system that can be solved for 
temperature and CO2 as a function of insolation (T(S) 
and (j>(S)). For simplicity we will take r = 1 in what 
follows, which is consistent with assuming a constant 
seafloor spreading rate and plate tectonics. 

3. EFFECT OF LAND FRACTION ON WEATHERING 

3.1. Simple Limits 

Throughout this paper we will be interested in the 
stability that the carbon cycle can provide to the cli- 
mate system. We will consider specifically the change 
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in surface temperature caused by a change in insolation 
which we will call the climate sensitivity to ex- 
ternal forcing. The lower S is, the more stable the 
climate is, the wider the habitable zone will be, and 
the longer a planet will stay in the habitable zone as 
its star evolves. We do not consider the many effects 
that could become important near habitab le zone lim- 
its, s u ch as water clouds on the inne r edge (jSelsis et al.l 
[2007t iKaltenegger fc Sasselovl [20Tll) and CQ 2 clouds 
on the outer e dge dForget fe Pierrehumbertj [T9971: 
ISelsis et al.l 120071: IWordsworth et al.l I201H iForget et al.l 
I2012t IWordsworth et al.l 120121) . but §f will give us a 
sense of how wide the habitable zone is in the broad 
range between where these effects could be important. 

We can calculate ^ by differentiating Equation (jlip 
with respect to S 



dT 
~dS 



LJ2 



4lui 



u>i(j) dS 



(12) 



The change in CO2 as the insolation changes (5^) de- 
fines the weathering feedback strength. In t he traditional 
pictu re of the silicate weathering feedback (|Walker et al.l 
119811 ). the CO2 decreases as the insolation increases 
(§3 < 0), stabilizing the climate. If ^ = there is no 
weathering feedback, and the temperature increases at a 
constant rate, = as the insolation increases. 

To calculate 5^, we need to consider weathering. We 
can generalize the weathering relations we described in 
Section [2] by rewriting Equation © as 



W = /(7)Wi*(^,T) + g^)W* {^T), 



(13) 



1 



— 1 

In 

in our stan- 



where /( 7 ) = 7, 3(7) = 1 - 7, W^T) 

a p Tf<j> b e T , and W*(^T) = Iff (^±f) ir 

dard model. The arguments in this section only require 
that /( 7 = 1) = .9(7 = 0) = 1 and /( 7 = 0) = .9(7 = 
1)=0. 

We will assume that the outgassing rate (r), and 
therefore the total weathering rate (r = W), does not 
depend directly on insolation (S). This is a reason- 
able assumption for surface temperatures in the hab- 
itable zone, although the chemical composition of out- 
gassed material does var y at much higher temperatures 
(jSchaefer fc Fegfevll2010D . Using this assumption, we can 
differentiate Equation (fT5|) with respect to S to find 



= /( 7 ) 



-9(l) 



ow; 



dcj> dS 
dW* dc 



dW t * dT 
dT 1)S 
d<f> dW* dT 

'as + 8t as 



(14) 



Equation ([14]) can be solved for which can be sub- 
stituted into Equation (|12[) to solve for IS. 

Before proceeding we consider two simple limits. First 
we consider a waterworld (7 = — > f = 0, g = 1) in 
which seafloor weathering does not depend on planetary 
surface temperature (W*(<j), T) = W* (</))), as implied by 
most seafloor weathering parameterizations (Section [2]). 



dW* dip 



In this case all that is left of Equation (|T4]) is —q^- 
If we assume that W* has a 4> dependence so that it can 



adjust to match outgassing, we are left with ^ = 0. We 
therefore see that the assumption that seafloor weather- 
ing does not depend on planetary surface temperature 
implies that there can be no weathering feedback on a 
waterworld. 
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Fig. 1. — Temperature and CO2 concentration as a function of 
insolation with modern Earth volcanic outgassing (r = 1) and var- 
ious land fractions: 7 = (dashed), 7 = 0.01 (light gray solid), 
7 = 0.1 (dark gray solid), and 7 = 0.99 (black solid). Both dimen- 
sionless (left) and dimensionful (right) values are given on axes. 
For reference, the time relative to zero age main sequenc e for a 
planet located 1 AU from a Sun-like star is plotted following lGougbl 
l|1981l ). Gray shaded regions of the plot represent conditions that 
may be considered uninhabitable by some definition. The habit- 
able zone upper temperature limit is the moist greenhou se, which 
we take to occur at 340 K (T = 5.3) IIKasting et al.H1993l) . and the 
lower temperature limi t is global glaciatio n, which we take to occur 
at 265 K (T = -2.2) dVoigt et al.l [20111). although it could occur 
10-20 K colder than this IIAbbot et aU20lil ; IYang et al.120121) . The 
lower limit placed on CO2 for habitability i n the figure, 10 -5 bar, 
llCaldeira fc Kasti ng 1992; Bloh et al.ll2fjQ5t) should be understood 
only as a potential limit on complex photosynthetic life, and not 
on other types of life. 

As a second limit, we allow the land fraction (7) to be 
non-zero but make seafloor weathering (W*) a constant 
independent of temperature and CO2. While this limit 
does not accurately describe our parameterizations in 
Section [5J we do expect seafloor weathering to be weaker 
than continental weathering an d to respond weakly to 
changes in CO2 (|Caldeiralll995l) . which makes this limit 
useful for gaining insight into the full model. Further- 
more, seafloor weathering may be limited by the avail- 
ability of basaltic crust for reaction, which would lead 
to constant seafloor weathering (Sections [2] and [7|). If 
seafloor weathering is set constant, then Equation (fT4| 
reduces to 



= m 



dW? d(j> dW t * dT 
dcf> OS + dT dS 



(15) 



0. Equation ((To"]) implies that the change in CO2 as inso- 



lation changes (§§) is independent of land fraction (7) 
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so that, by Equation (jT2J) , the climate sensitivity to ex- 
ternal forcing must also be independent of land 
fraction. This important result will lend insight to the 
full model, in which we relax the assumption of constant 
scafloor weathering, where we find that ^ is insensitive 
to land fraction (7), even for very small land fraction. 
Note that as long as f(<y) is monotonically increasing 
and all other parameters are held fixed, the continental 
weathering rate decreases as the land fraction decreases 
(Equation (fT3| ). This would lead to increases in temper- 
ature and CO2 (T and </>) as land fraction (7) decreases, 
but they adjust precisely to keep ^ constant in the limit 
of constant seafloor weathering. 

3.2. Using the Full Model 

We now consider the full model described in Section [5] 
rather than just the simple limits considered in Section 
13.11 This means, for example, that seafloor weathering 
is given by Equation (J7|) rather than held constant. In 
the full model there is a striking difference between the 
weathering behavior of a waterworld and a planet with 
even a little land (Figure [TJ). Varying the land fraction 
(7) from 0.99 to 0.01 causes only a moderate warming 
and has a very small effect on the climate sensitivity to 
external forcing (§^)- The functioning of the silicate 
weathering feedback in the partially ocean-covered plan- 
ets can be seen from the decrease in CO2 (<t>) with insola- 
tion (S; Figured]). Decreasing 7 from 0.99 to 0.01 would 
therefore shift the habitable zone outward in insolation 
(S) slightly, but would not significantly alter its width. 

Setting the land fraction to zero (7 = 0), however, 
causes a huge shift in planetary weathering behavior. 
With 7 = 0, the CO2 builds up to a high enough value 
that scafloor weathering can accommodate outgassing, 
and stays at this high value regardless of insolation (S 1 ), 
so there is no weathering-climate feedback. This leads 
to a large increase in surface temperature, CO2, and the 
climate sensitivity to external forcing (§^) for the water- 
world relative to the partially ocean-covered planets, so 
the waterworld crosses the moist greenhouse threshold 
much sooner (Figure [JJ . Therefore a planet with even 
a small amount of continent should enjoy a much wider 
habitable zone, and longer tenure in it, than a water- 
world. 

None of the curves with 7 > in Figure [1] cross the 
moist greenhouse threshold for S < 1.4 due to the sim- 
plifications in radiative transfer we have made, which is 
unrealistic when com pared to more detailed calculations 
(jKasting et al.1 fl993T k however, it is clear that a planet 
with a lower land fraction will tend to be warmer and 
experience a moist greenhouse at a lower insolation (S). 
On the other hand, a lower land fraction (7) leads to a 
higher CO2 (</>), so a planet with a lower land fraction (7) 
will cross the CO 2 (<fi) threshold for C4 photosynthesis 
(jBloh et al.ll2005l ) at higher insolation (S). 

We can solve for climate sensitivity to external forcing 
(f^) explicitly as follows 

dT (l-a\( & + e(T,</),7) \ 

dS l^rj [^(l + T ^) + b + Q(T,cj >ll )) ' 

(16) 



b=0.1 




Fig. 2. — Ratio of the climate sensitivity to external forcing (^) 
to the value it takes if there is no weathering feedback = 

4 ~ Q ), as a function of the insolation (S) and the land fraction 

(7) with 6=0.1 (upper panel) and 6=0.5 (lower panel). 6 is the 
exponent of atmospheric CO2 in the continental weathering rate 
parameterization (Equation |T5J| )■ Our standard value of b is 0.5, 
and we consider a lower value because the effect of atmospheric 
CO2 concentration on continental weatherin g may be reduced if 
land plants are present (fPicrrchumbcrt 20ToT) . 

where 

/ (j) n - b e - T \ 

X U + ")(l + °pT)°J ' (17) 

Note that the factor multiplying = (§^) in Equa- 
tion (|16|) is always less than or equal to one, which re- 
duces the climate sensitivity to external forcing {^) be- 
low (f^) and is a statement of the silicate weather- 
ing feedback. Also, if /3q = 1 (no scafloor weathering), 
9 = and ^ has no dependence on the land frac- 
tion (7), as we expected from Section [XT] Additionally, 
lim 7i/3o ^ = oc and lim e ->oo §§ = (ff) 7 which re- 
covers the waterworld limit. 

We plot the climate sensitivity to external forcing ( ^ , 
Equation (fT6"]) ) as a function of insolation (S) and land 
fraction (7) for different values of b in Figure [2] It is 
useful to consider a lower value of b because the effect 
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of CO2 (0) on continental weather ing may be greatly re - 
duced if land plants are present (|Pierrehumbertl [2(3101 ) . 
as we will discuss further in Section [4] The climate sen- 
sitivity to external forcing (§f ) increases for small land 
fraction (7), but insolation (S) has a much stronger effect 
on ^ than land fraction (7) does. The climate sensitiv- 
ity to external forcing (|f) is highest at low insolation 
(S) and land fraction (7) and decreases as S and 7 in- 
crease. lim 7 ^i G = and 9 — > as S becomes large 
because T becomes large (Equation (jTT)0 . This allows 
the approximation that if either the planet is more than 
«10% land covered and/or the insolation is as large as 
or greater than the modern Earth value, ^ approaches 
the limit 



dT 
~dS 



&(1 



dT 
dS 

dT 



l+a„T 



) + & 



(18) 



Therefore, despite the many parameters in our model 
(Section [2]), the climate sensitivity to external forcing 
(§|, Equation (|TB)) ) can be reduced to a function of only 
a few relevant variables over most of the land fraction- 
insolation plane (Equation (|18|l ). In particular, if we ne- 
glect the (usually small) term jt^~t > * ne om y weathering 
parameter present in Equation (fT8"|) is b, which controls 
the strength of the scaling of land weathering with CO2 
(other remaining variables are from the climate model). 
The weaker land weathering scales with CO2 (lower b), 
the less T can change while maintaining a balance be- 
tween outgassing and weathering, and the lower is. 
In the context of this paper, what is important about 
Equation (fT8| is that for a broad range of insolation, 



the climate sensitivity to external forcing (^) is near 
the limit given by Equation (|18[) and is therefore nearly 
independent of land fraction (7). This explains the in- 
sensitivity of the climate sensitivity to external forcing 
(§f ) to land fraction in this model. 

4. SENSITIVITY TO MODEL PARAMETERS 

In this Section we will evaluate the robustness of our re- 
sult from Section [3] that the weathering behavior, specif- 
ically the climate sensitivity to external forcing (S), 
of a planet is relatively insensitive to surface land frac- 
tion (7). We will compare ^ for a small land frac- 
tion (7 = 0.01) with that for an Earth-like land frac- 
tion (7 = 0.3) while we vary model parameters. For 
simplicity we consider reference insolation, outgassing, 
and weathering values (S = T = W = 1) and define 

A = (H) |t=0,^=1,S=1,7=0.01 - (§§) |t=0,^=1,ST=1,7=0.3 

as a measure of the change in |f for a large change in 
land fraction. 

We first consider the effect on our results if current es- 
timates of the fraction of total weathering accomplished 
on continents on modern Earth (/3 ) were incorrect. A 
lower /?o would increase |f (Equation (fTBjl). and this ef- 
fect is exaggerated when the land fraction (7) is smaller, 
increasing A (Figure [3]). For /3o less than approximately 
0.5, A starts to become relatively large. Given current 



estim ates of 0.67 (|Caldeiralll995h and 0.78 (|Le Hir et al.l 

120081) for /3o, seafloor weathering on modern Earth would 
have to be 2-3 times larger than it is thought to be to 
affect our main conclusions. 

Next we consider changes in the seafloor weathering 
scheme. Ionic conce ntration in sea water may either in- 
crease or decrease n (jCal deira 19951) , the exponent of 4> in 
the seafloor weathering parameterization (Equation ([7])). 
The constant m in Equa tion (|6l), howev er, is typically 
taken to be less than 1 (jCaldeiral 119951 ). which would 
correspond to n = |. Such a change would lead to a 
moderate increase in A (Figure [3]), although a larger in- 
crease in n would significantly increase A. It is less clear 
whether a vastly different value of v could be appropriate, 
but only an extremely large decrease in v would signif- 
icantly increase A (Figure [3]). The large effect that /3o 
and n would have on A if appropriate values were out- 
side current ranges highlights the importance of more ex- 
perimental and field work to better understand seafloor 
weathering. 

The exponent of precipitation in Equation ((3]), a, and 
the constant defining the increase in precipitation with 
temperature (Equation (O), a p , have a small effect on §^ 
(Figure ■ This is because the exponential temperature 
dependence of the Arrhenius component of the conti- 
nental weathering parameterization (Equation (j3])) dom- 
inates changes in continental weathering due to temper- 
ature. Moreover, changes in a and a p affect ^ similarly 
at low and high land fraction (7), leading to only very 
small changes in A when these parameters are changed 
(Figure [21). 

We next consider changes in 6, the exponent of CO2 in 
the continental weathering parameterization. If seafloor 
weathering is constant , setting b = yields a perfect 
thermostat (^ = 0), since weathering is a function of 
T alone in this case (Equation ([8])). As b is decreased 
in the full model, |f decreases, with larger changes for 
larger land fraction (7), since the model behaves more 
like the constant seafloor weathering case in this limit 
(Figure 13]) ■ For very small b, A becomes somewhat large, 
although this effect is greatly diminished for insolation 
larger than the reference value S > 1 (Figure [2]) , when 
the temperature becomes large enough to make small 
(Section [3]). 

Finally, we consider the terms defining the expansion 
of infrared emission to space (f2, outgoing longwave ra- 
diation) as a function of temperature (wi = |p) and the 
logarithm of CO2 (w 2 = ~ diog(0) )- ^ e on ^ kept ^ ne nrs t 
order terms in these expansions, and we made specific as- 
sumptions to calculate them (Section so it is worth 
considering whether variations in lo\ and uji affect our re- 
sults. Decreasing uj\ has the simplistic effect of increasing 
the climate sensitivity to external forcing with no weath- 
ering feedback (§f ) = "3^7 • We focus on weathering 
feedbacks in Figure |3]by normalizing the climate sensitiv- 
ity to external forcing by (§f ) - and u>2 enter the re- 
lation for the weathering feedback in Equation (fTS)) (the 
part multiplying (§f ) on the right hand side) through 
the ratio — , which is a measure of the climate sensitiv- 
ity to changes in CO2. Increasing the climate sensitivity 




Fig. 3. — Ratio of the climate sensitivity to external forcing (S?) to the value it takes if there is no weathering feedback (( §S 



l — a \ 
4 Wl > 

at both low land fraction (7 = 0.01, thick black solid line) and moderate land fraction (7 = 0.3, thick gray solid line). ^ is displayed as 
a function of /3o, the fraction of total weathering accomplished on continents on Earth today; n, the exponent of atmospheric CO2 (</>) in 
the seafloor weathering parameterization (Equation 10); v, the constant in the seafloor weathering parameterization; a, the exponent of 
precipitation in the continental weathering parameterization (Equation 0); ce p , the constant defining the increase in precipitation with 
temperature (Equation 10); b, the exponent of atmospheric CO2 (</>) in the continental weathering parameterization (Equation 0); uii, 
the first term in the expansion of the outgoing longwave radiation (f2) with temperature (jjp; Equation flo]) ); and Cd2, the first term in 



the expansion of outgoing longwave radiation with the logarithm of CO2 (- 
standard values of parameters. 



£>iog{« 



Equation mOH V Circles and the thin black line signify 



to changes in CO2 (— ) increases the warming effect of 
a given increase in CO2, increasing the weathering feed- 
back strength, and decreasing the climate sensitivity to 
external forcing (Figure [3]) . This effect shows little de- 
pendence on land fraction, so large changes in u)\ and L02 
do not alter our main results (Figure [3]) . 

We note that T u , the exponential temperature scale 
for the direct temperature dependence in Equation ([2]), 
comes into the nondimcnsional relations through a p = 

a p T u and u>i = wiif-- Since reasonable variations in a p 
and oj\ do not alter our main conclusions, we can con- 
clude that reasonable changes in T u would not. Indeed, 
we find that the climate sensitivity to external forcing 
depends minimally on land fraction (7) when we vary T u 
between 5 and 25 K (not shown). 

5. WATERWORLD SELF-ARREST 

Due to the lack of a climate-weathering feedback in 
a waterworld, we have argued that a waterworld should 
experience a moist greenhouse at much lower insolation 
(S) and have a smaller habitable zone than a planet with 
even a very small land fraction. We will now discuss the 
possibility that a waterworld could stop a moist green- 
house by engaging silicate weathering when enough ocean 
is lost that continent is exposed, engaging weathering 
and reducing temperature (T), and thereby preventing 
complete water loss. We will refer to this possibility as 
a "waterworld self-arrest" (see Figure [4] for a schematic 
diagram of the idea). The waterworld self-arrest is sim- 



ilar to the idea, discussed bv lAbe et al.l (|2011f ). that a 
planet with partial ocean coverage could experience a 
moist greenhouse and become a habitable land planet. 

If a waterworld experiences a moist greenhouse, there 
will be some amount of time before any continent is ex- 
posed and continental weathering is possible. We call 
this timescale tq (Figure [5]) , and we will not use it in the 
present argument. We call the time after the first expo- 
sure of continent until all remaining water is lost through 
hydrogen escape to space r rn . To determine whether a 
waterworld self-arrest is possible, we must compare r m 
to the timescale for CO2 to be reduced by continental 
silicate weathering to a level such that moist greenhouse 
conditions cease to exist, which we call t w . If t w < t to , 
then a waterworld self-arrest is possible. 

If we take the hydrogen atom mixing ratio in the 
stratosphere to be unity as a limiting case, the diffusion- 
limited rate of loss of hydrogen atoms to s pace is 2 x 10 13 
H atoms cm" - 2 s- 1 (jKasting et al.lll993j ). Earth would 
be essentially a waterworld if the ocean were a few times 
more deep, which would correspond to « 10 29 H atoms 
cm~ 2 , so that r m is of order 100 Myr. Since we seek only 
an order of magnitude here, we do not consider the effect 
of different topography on r m , but note that larger plan- 
ets will tend to have less topography, which would reduce 
T m somewhat. We note also that en ergetic limitation on 
hydrogen escap e could increase r m (| Watson et al.lll98lT ; 
I Abe et alJ[201lh . 

t w is much more difficult to estimate than r TO . The 
closest known analog to a moist greenhouse in post- 



10 



Abbot, Cowan, & Ciesla 




340 



310^ 




280 



250 



Fig. 4. — Schematic diagram showing how the waterworld from 
Figure[T]could experience a waterworld self-arrest as the insolation 
(5) increases and end up as an Earth-like planet with partial ocean 
coverage Oy = 0.3 in this case) in the habitable zone. Axes are as 
in Figure [T] 



Waterworld Self-Arrest Timescales 

water level 



water loss 



initially 
submerged 
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T m : time to complete 
water loss after 
continent exposed 



Fig. 5. — Schematic diagram of the timescales involved in a wa- 
terworld self-arrest, to is the time the planet exists in a moist 
greenhouse state before submerged continent is exposed. r m is the 
time to complete planetary water loss after the first exposure of 
continent. In order to determine whether a waterworld self-arrest 
is possible. r m must be compared with t w , the timescale for CO2 
to be sufficiently reduced by silicate weathering that moist green- 
house conditions cease to exist. 

Hadean Earth is likely the aftermath of the Snowball 
Earth events, which occurred 600-700 million years ago. 
The dominant explanation for Snow ball Earth events 
(|Kirschvin]dll99aiHoffman et al.lll998l ) is that some cool- 
ing factor forced Earth into a global glaciation, which is 
an alternative stable state of climate because of the high 
reflectivity of snow and ice. Because it would be very cold 
during a Snowball and because any precipitation would 
fall as snow, rather than rain, the continental weather- 
ing rate would be near zero. As we expect volcanic CO2 
outgassing to continue unabated, the atmospheric CO2 
would increase until it reached up to a thousand times 
its present value (0 = O(10 3 )), at which point it would 
be warm enough to melt the ice and end the Snowball. 
In the immediate aftermath of a Snowball Earth event 



the planetary albedo would be very low and the CO2 
concentration would be very high, which would lead to 
surface temperatures in the range of a moist greenhouse. 
Extremely high weathering would then reduce the at- 
mospheric CO2 until Earth's s urface temperatures wer e 
comparable to modern values (jHiggins fc Schragll2003| ). 
If post-Snowball Earth was truly in a moist greenhouse 
state, complete water loss must not have occurred, which 
would imply that t w < r TO and a waterworld self-arrest 
is possible. Even if post-Snowball Earth was not a true 
moist greenhouse, the drawdown of CO2 must have been 
immense, and geochronolog ical evidence indicat es that 
this took less than 4 Myr (jCondon et al.ll2005[ ). This 
further motivates the idea that t w could be small com- 
pared to T m . 

The CO2 required to end a Snowball Earth, however, 
may have been 1-2 orders of magnitude lower than that 
quoted above (jAbbot fc Pierrehumbertl 12010( 1 . in which 
case the Snowball aftermath would be a weaker analog 
for a moist greenhouse and further consideration is re- 
quired. Geochemical modeling under conditions similar 
to a moist greenhouse yield weathering rates of 0(1O 13 ) 
kg C yr- 1 (IHiggins fc Schrad l200l and 0(1O 12 ) kg C 
vr" 1 (|Le Hir et all l2009h . although [MnisetalJ IpOll 
suggest that ion transport constraints coul d limit the 
weathering rate to « 5 x 10 11 kg C yr -1 ([Mills et al.l 
120111 1 . The waterworld C0 2 is 1.65 x 10 4 times the 
present value (Figure [T]) in our model with standard pa- 
rameters, which corresponds to O(10 19 kg C). This yields 
an estimate of t w — 0(1-10 Myr), which further suggests 
that t w <C Tjyi . 

Although we have not performed a full analysis of the 
kinetic (non-equilibrium) effects, the order-of-magnitude 
analysis we have done indicates that a habitable zone wa- 
terworld could stop a moist greenhouse through weath- 
ering and become a habitable partially ocean-covered 
planet. We note that this process would not occur if 
the initial water complement of the planet is so large 
that continent is not exposed even after billions of years 
in the moist greenhouse state (t > O(10 Gyr)). Addi- 
tionally, we have neglected the possibility of water ex- 
change between the mantle and the surface, which is 
poorly constrained on Earth. If, for some reason, water 
were continuously released from the planetary interior 
to replace water lost through hydrogen escape, then the 
conclusions reached here would be invalid. This would 
require significantly higher mantle-surface fluxes than on 
present Earth, which are constrained by cerium mea- 
surements at mid-ocean ridges to be small enough to 
requi re ~10 Gyr to replace Earth's s urface water reser- 
voir (jHirschmann fc Kohlstedtl 12012( 1 . Finally, we note 
that since we have considered a moist greenhouse caused 
by high atmospheric optical thickness rather than high 
insolation, the waterworld would not be likely to suffer 
a full runaway greenhouse during the process described 
here (|Nakaiima et alj|1992l) . 

6. DISCUSSION OF OBSERVABILITY 

The theory developed here suggests that if an Earth- 
size planet is detected near the habitable zone it should 
be considered a good candidate to have a wide habitable 
zone and a long-term stable climate as long as it has some 
exposed land. The first observational point to consider, 
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therefore, is whether the surface land fraction could be 
measured on an Earth-like planet. 

Using measurements of the rotational variations of 
disk-integrated reflected visible light from a future mis- 
sion such as TPF-C it will be possible to determine the 
surface land f raction of an Earth- like planet orbiting a 
Sun-like star. iCowan et al.l (|2009f ) used observations of 
Earth from the Deep Impact spacecraft to demonstrate 
that this is possible with multi-band high-cadence pho- 
tometry. The multi-band observations are critical be- 
cause the albedo of land and ocean have very differ- 
ent wavelength dependencies, while the integration times 
must be sufficiently short to resolve the planet's rota- 
tional variability. Clouds tend to obscure surf ace features 
and make such a determination more difficult dFord et al.l 
12001 iKaltenegger et al.ll2007t iPalle et alJl2008f ). and re- 
gions of a planet that are consistently cl oudy are im- 
pervi ou s to this form of r emote sensin g (jCowan et al.1 
[20ll . iKawahara fc"Fuiiil (|2010l I20H have extended 
this method to include seasonal changes in illumina- 
tion to generate rough two-dimensional maps of other- 
wise unresolved planets. Provided that a habitable zone 
planet is not much cloudier than Earth, its land fraction 
should therefore be observationally accessible with suffi- 
cient signal-to-noise ratio, and the present study can be 
used to help assess the probability of its habitable zone 
width and long-term habitability. 

The presence of exposed land could be detected on 
an cxoplanct via an infrared emission spectrum, pro- 
vided the overlying atmos phere is not too much of a 
hindrance (|Hu et al.l I2012T ) . At thermal wavelengths, 
the land fraction of a planet could in principle be es- 
tima ted based on its thermal inertia, but in idealized 
tests lGaidos fe Williams! (|2004l ) found that the obliquity- 
inertia degeneracy impedes quantitative estimates of a 
planet's heat capacity. It is therefore unlikely that a 
thermal mission, on its own, could constrain the land 
fraction of an exoplanet. 

Direct confirmation or falsification of the theory de- 
veloped here will be more difficult, requiring not only 
an estimate of land fraction, but also atmospheric CO2 
abundance and/or surface temperature. The difference 
between the temperature and CO2 of a waterworld and a 
partially ocean-covered planet could be quite large, par- 
ticularly near the inner (hotter) edge of the habitable 
zone (Figure [T]). This entails measuring, in exoplan- 
ets, discrepancies in carbon dioxide abundance of roughly 
four orders of magnitude, and surface temperatures dif- 
ferences of <~ 70 K. 

Detecting carbon dioxide in an cxoplanct atmo- 
sphere is relatively easy: the James Webb Space Tele- 
scope will determine whether CO2 is present in the 
atmos pheres of transiting super-Earths orbiting M- 
stars (iDeming et all 120091: IKaltenegger fc Traubl 120091 : 



iBelu et al.1 1201 ll iRauer et alll2011[) . The challenge is in 
estimating the amount of CO2 in an atmosphere. In or- 
der to measure the carbon dioxide abundance of an atmo- 
sphere, the atmospheric mass must be estimated. This 
could be achieved for M-Dwarf habitable zone planets 
with tran sit spectroscopy covering the Rayleigh scatter- 
ing slope (jBenneke fc Sc ager 2013). The unknown radius 
of directly-imaged (longer-period) cxoplancts may hinder 
such estimates with a TPF-C type mission. At thermal 
wavelengths, the planetary radius will be measurable and 



the pressure broadening of absorption features may offer 
an estimate of surface pressure and hence atmospheric 
mass (e.g., iDes Marais et al.l 120021 : iMeadows fc Seagerl 
[20101) . 

Using a thermal mission like TPF-I, it may be possible 
to determine the surface temperature of an Earth ana- 
log via its brightness temperature in molecular opacity 
windows, at l east for a cloud-free atmosphere with few 
PR absorbers (IDes Marais et alj [20021 : IKaltenegger et all 
120101 : iSelsis et all 120 111 ). While low-lying clouds do not 
significantly degrade such estimates, high cloud coverage 
can lead to surface temperature estimates at least 35 K 
too low (jKaltenegger et al J 120071: iKitzmann et alll2011ft . 
Continuum opacity from water vapor will also lead to 
lower brightness temperatures, even in opacity windows. 
This could lead to under-estimating the surface temper- 
ature of planets with high vapor content, precisely the 
moist greenhouse planets expected to have high surface 
temperatures. Moreover, Venus dramatically demon- 
strates that brightness temperatures, even in supposed 
windows, are not necessarily a good probe of the surface 
temperature. 

Given the similar temperature and CO2 among par- 
tially ocean-covered planets of varying land fraction (Fig- 
ure [T|) , it is unlikely that any mission in the foreseeable 
future would be able to distinguish between them. Ad- 
ditionally, factors that we have neglected here introduce 
uncertainties that make the temperature and CO2 abun- 
dance among partially ocean-covered planets effectively 
indistinguishable . 

To summarize, measurements from a mission like TPF- 
C would provide a surface land fraction estimate for plan- 
ets not much cloudier than Earth. TPF-C could pro- 
vide an order-of-magnitude atmospheric CO2 estimate 
if the radius were known, which would be possible if 
(1) the planet were transiting, (2) its radius were esti- 
mated from TPF-I observations, or (3) its radius were 
estimated based on a mass measurement. TPF-I could 
provide a surface temperature estimate that would likely 
be good to O(10 K) for an Earth analog, but possibly 
much worse for a cloudier and/or more humid atmo- 
sphere. It would therefore be possible to directly test 
our predictions for differences in weathering behavior of 
waterworlds and partially ocean-covered planets if data 
from missions like both TPF-C and TPF-I were available, 
and sufficient numbers of nearby habitable zone planets 
exist. Prospects are limited, but testability is still possi- 
ble in special cases (such as a transiting planet), if only 
one type of data is available. 

A final observational point is that in Section [5] we cal- 
culated that a waterworld self-arrest, through which a 
fully ocean planet would be transformed into a partial 
ocean planet, is possible. Waterworlds, however, would 
still exist and their detection would not falsify the water- 
world self-arrest idea. Waterworlds coul d have orders of 
magnitude more wa ter than Earth (e.g.. iRavmond et all 
120071: IFu et alll2010[ ). and could therefore exist in a moist 
greenhouse state for many billions of years. Additionally, 
a waterworld in the outer regions of the habitable zone 
would not be in a moist greenhouse state (Figure H|) . 

7. GENERAL DISCUSSION 

We note again that we make numerous approximations 
in our model, so that it is useful for qualitative under- 
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standing rather than quantitative predictions. Our lin- 
earized climate model, for example, is a gross simplifi- 
cation of the real climate system, which is probably re- 
sponsible for the delayed onset of the moist greenhouse 
in our model as the insolation is inc reased (FigureJ H) rel- 
ative to more detailed calculations (jKasting et al.lll993h . 
Our approach, however, is appropriate for a general 
study of potential terrestrial exoplanets, whose detailed 
characteristics are not known, and our main conclusions 
should be robust to more detailed modeling. For exam- 
ple, we point out that no weathering feedback is possi- 
ble = (f^) ) for a waterworld as long as seafloor 
weathering is independent of surface temperature (T), 
regardless of the specific form of the seafloor weathering 
parameterization. 

Nonetheless, our parameterization of seafloor weather- 
ing is probably the most uncertain aspect of this work. 
For example, we have neglected the potential effect of 
a different ocean chemical composition on an cxoplanet, 
which could lead to different base weathering rates or 
a different scaling of seafloor weathering with CO2 (</>)■ 
Furthermore, in order to derive a scaling relationship be- 
tween [H + ] and 0, we have approximated ocean chem- 
istry as a carbonate ion system with changes in alkalinity 
determined by changes in Ca 2+ . Real oceans have other 
ions present, which can alter the charge balance. These 
approximations, however, have been useful for model- 
ing carbon cycling in deep time Earth problems (e.g., 
IHiggins fc Schral200l and our knowledge of typical ion 
concentrations in terrestrial exoplanets is unconstrained. 
Finally, we have assumed that surface temperature has 
a negligible effect on seafloor weathering rate, which we 
will now consider in more detail. 

As there is a temperature depen dence of the reac- 
tions relevant for seafloor weathering (jBradv fc Gislasonl 
119971 ). we have effectively assumed that surface temper- 
ature does not strongly affect the temperature at which 
seafloor weathering occurs. If a dependence of seafloor 
weathering on surface temperature is allowed, a climate- 
seafloor weathering feedback is possible, which should 
tend to reduce the effect of changing land fraction. In 
terms of the effect of land fraction on weathering behav- 
ior, we have therefore considered a conservative case in 
assuming no surface temperature dependence of seafloor 
weathering, yet we still found that land fraction has only 
a small effect on weathering behavior for partially ocean- 
covered planets. 

A strong temperature dependence of seafloor weather- 
ing would affect waterworld behavior at the inner edge of 
the habitable zone. We can see this by multiplying the 
seafloor weathering parameterization (Equation 0) by 
an exponential temperature dependence of the form e r,T , 
where r\ = 4^ , where d is the change in seafloor weather- 
ing reaction temperature per unit change in surface tem- 
perature, and T w w 15K is the exponential scale of in- 
creas e in seafloor weathering reaction rate with tempera- 
ture (jBradv fc Gislasonlll997[) . The largest value r\ could 
take would be if current reactions occur near the bot- 
tom water temperature («0°C) and each unit increase 
in surface temperature creates a unit increase in bottom 
water temperature (d=l, -q — |). Our standard assump- 
tion is that seafloor weathering reactions occur at a tem- 
perature of p»20-40°C regardless of surface temperature, 



which is equivalent to d=j]=0. Using rj = | increases the 
insolation required to cause a waterworld moist green- 
house from 5«1.2 to S~2.1 if v~Q and prevents a water- 
world moist greenhouse from ever happening if ^=12.6. 
A strong temperature dependence of seafloor weathering 
would therefore make it much harder for a waterworld 
to experience a moist greenhouse, although we do not 
consider the limiting case considered in this paragraph 
likely. 

As noted in Section [21 we assume that seafloor spread- 
ing produces a sufficient quantity of weatherable material 
and a sufficient amount of water can circulate through it 
that the seafloor weathering rate can increase arbitrar- 
ily as the ocean pH decreases. Since we assume scafloor 
weathering accounts for 25% of weathering on modern 
Earth (/3q = 0.75), if the CO2 outgassing rate takes its 
modern value (r = 1), then the seafloor weathering rate 
must be able to increase by a factor of four if it alone is to 
balance outgassing, as would be required to establish a 
CO2 steady state on a waterworld. If there is not enough 
weatherable material for this to be possible, for example 
due to weak circulation of seawater through weatherable 
material, then the waterworld CO2 will not be limited 
to the value calculated in Section [3] (Figure [TJ), but will 
instead steadily increase due to an unbalanced source 
and sink. This would push a waterworld more quickly 
to a moist greenhouse, and potentially a waterworld self- 
arrest, than one would otherwise expect. A related issue 
is that if the land fraction is very small, there may be 
physical limits on continental weathering not included in 
Equation ([3]). This would lead to higher CO2 concen- 
trations and temperatures for small land fractions than 
those found in Section [3] The smallest land fraction we 
explicitly considered, however, was 7 = 0.01, which cor- 
responds to the combined area of Greenland and Mexico 
on Earth. Given continuous resurfacing it is not clear 
that physical limitations on weathering would be in play 
on such a large land area. 

We have focused on the effect on land fraction on 
weathering behavior relevant for the broad interior of the 
habitable zone and the evolution of a planet through the 
habitable zone over time. The effect of land fraction on 
the outer (cold) edge of the habitable zone is of interest 
as well. Based on increased Rayleigh scattering or CO2 
condensation at the surface, the outer ed ge of the hab- 
itable zone occurs when C O2~5-10 bar (jKasting et al.l 
119931 ). iLe Hir et aH (|2008| ). however, found that seafloor 
weathering should limit the maximum CO2 concentra- 
tio n on Earth to a few t enths of a bar. As pointed out 
by iPierrehumbertl (|2011[ ). if this is the case, then the 
CO2 condensation limit could not be reached and seafloor 
weathering would shrink the habitable zone. Using our 
seafloor weathering parameterization, we can derive an 
upper limit on the atmospheric CO2 concentration with 

Earth's land fraction (7 = 70): <j> = ( j^zrj^ — v> ) ■ Us- 
ing our standard parameters, we find <j> = 7.3 x 10 4 , or 
«20 bar. cf> depends strongly, however, on model param- 
eters (Figure El), particu larly n. <fi = 10 3 , roughly the 
value iLe" Hir et all (j2008|) found, would be produced by 
our model by either a change in n from 0.33 to 0.54 or 
a change in v from 12.6 to 2.0. Therefore, in contrast 
to our main results (Section 0]), the weathering behavior 
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relevant for the outer edge of the habitable zone is highly 
sensitive to the seafloor weathering parameterization. A 
convincing study of the effect of seafloor weathering on 
the outer edge of the habitable zone will have to wait 
until seafloor weathering is better understood. 



log 10 (CO™ x ) [bar] 




0.2 0.4 0.6 0.8 1 



n 

Fig. 6. — Maximum atmospheric partial pressure of CO2 (dimen- 
sionful version of <f>) allowed by seafloor weathering, as a function 
of n and v, parameters in the seafloor weathering scheme (Equa- 
tion 0). The "+" represents standard parameter values. 

We have found that relatively ineffective seafloor 
weathering will tend to lead to much warmer temper- 
atures in a waterworld than a p artial ocean plan e t. Thi s 
might appear to contrast with iSleep k Zahnld (|2001| ). 
who calculated very cold temperatures for Earth during 
the Hadcan, when Earth had a much lower land fraction 
and was further out in the habit able zone. The p r imary 
reason for the difference is that ISleep k Zahnlei (|2001l ) 
include a significant sink of carbon due to weathering 
of impact ejecta during the Hadean, which significantly 
reduces the equilibrated CO2 and temperature. They 
also assume that seafloor weathering is much more ac- 
tive due to increased internal heat flux early in Earth's 
history. We do not include either of these effects be- 
cause our effort is aimed at general consideration of the 
problem in the context of exoplanets, but we acknowl- 
edge that either could limit the CO2; and therefore the 
proba bility of a moist greenh ouse, on a waterworld. Fi- 
nally, ISleep k Zahnlei ([200 ll ) also find a strong depen- 
dence of equilibrated CO2 on the strength of the CO2 
dependence of seafloor weathering, as we do (Figure [5]), 
early in Earth's history when the CO2 may have been 
very high. This again emphasizes that the details of the 
seafloor weathering scheme are critical for its effect on 
the outer edge of the habitable zone, or equivalently the 
climate of early Earth. 

In our analysis we set the planetary albedo to a con- 
stant for simplicity. Since the albedo of land (0.05-0.3, 
IPeixoto ~k Oort 19 92) tends to be higher than the albedo 
of ocean (0.05-oT TPeixoto k Oortj[l992T) , one would ex- 
pect the planetary albedo to increase as the land frac- 
tion increases. This would decrease the temperature (T, 
Equation (|TT|) ) and the climate sensitivity to external 



forcing (^, Equation ([TB")) ). both of which we found 
would already be lower due to more continental weath- 
ering (Section [3]). This effect would likely not be very 
large, however, as detailed analysis of satellite measure- 
ments indicates that because of extensive absorption and 
reflection by clouds and gas in the atmosphere the sur- 
face contri butes only 0(10%) to the global mean albedo 
on Earth (|Donohoe k Battistill201lD . The influence of 
surface albedo on planetary albedo would be further re- 
duced around M stars, which emit at longer wavelengths, 
due to increased atm ospheric scattering and absorption 
(jKasting et al.llT993D . Another surface albedo effect we 
have neglected is the change in albedo in colder climates 
due to the formation of ice and snow at the surface. This 
would lead to a reduction in the range of S for which a 
planet is habitable in Figure [TJ Again, however, this 
effect would be reduced aroun d M stars due to grea tly 
reduced ice and snow albedos (jJoshi k H abcrlc 2 01lT ). 

Our model neglects spatial variation in climate and 
weathering, and therefore cannot address spatially het- 
erogeneous effect s such as the feedbacks described by 
Kite et all (|2011l ) . An assumption of spatial homogene- 
ity, however, is a reasonable first approximation even on 
planets with highly asymmetric insolation distributions, 
such as tidally locked planets in the habitable zone of M 
stars, if there is an ocean or atmosphere si milar in size 
to those on Earth to redistribu te heat (e.g.. JJoshi et all 
[1991 iMerlis k Schneide3l2010t lEdson et al]l20ll ~ 

Current understanding is that continents were not 
present at Earth's formation, but started to form 
soon afterward. For example, zircon evidence indi- 
cates that there was continental crust 4. 4 Gyr ago, 
or wi thin wlOO Myr of Earth's formation ([Wilde et al.l 
120011 ). Available evidence cannot constrain whether 
the continents reached their present size quickly (within 
«500 Myr) or have grow n slowly over Earth's h istory and 
continue to grow today (jEriksson et al.l [20061 ). We have 
neglected such processes, but this is unlikely to mod- 
ify our main conclusions because we found that the cli- 
mate sensitivity to external forcing (§^) does not de- 
pend strongly on land fraction. Our work suggests that 
the order of magnitude of water delivery, which deter- 
mines whether a planet is dry, partially ocean-covered, 
or a waterworld, is much more important for determining 
weathering behavior than evolution in the size of conti- 
nents through a planet's lifetime. Finally, it is worth 
speculating that Earth could have experienced a water- 
world self-arrest early in its history, although much more 
thought is needed to determine what the geochemical 
signals left behind by such an event would be. 

If a waterworld has a thick enough ocean, the overbur- 
den pressure could significantly suppress volcanic out- 
gassing. For example, CO2 degassing could be com- 
pletely s uppressed on an Earth-like planet with 100 km 
of ocean ([Kite et al.ll2009D . This uncertainty is commen- 
surate with other volcanic outgassing uncertainties, such 
as CO2 content in the magma, for partial ocean planets 
and waterworlds that barely cover their continents with 
ocean. Overburden outgassing suppression could lead to 
lower CO2 levels and higher chance of a Snowball state 
than predicted here for waterworlds with a thick ocean 
layer. 

In our discussion of a waterworld self-arrest (Section[5]) 
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we neglected effects of changing planetary size. Super- 
Earths would have higher gravity and therefore reduced 
orography. This will reduce the amount of time it takes 
from first exposure of land to complete water loss (r m ), 
which will make a watcrworld self-arrest less likely. On 
the other hand, the higher gravity of a super-Earth will 
also make gas escape from the atmosphere to space more 
difficult, which would reduce r m . Full consideration of 
a waterworld self-arrest on a super-Earth is beyond the 
scope of this paper. 

8. CONCLUSIONS 

We have shown that using standard weathering pa- 
rameterizations, the weathering behavior of a partially 
ocean-covered Earth-like planet does not depend strongly 
on land fraction, as long as the land fraction is greater 
than «0.01. Consequently, planets with some continent 
and some ocean should have a habitable zone of simi- 
lar width. This is a powerful result because it indicates 
that previous habitable zone theory developed assuming 
Earth-like land fraction and weathering behavior should 
be broadly applicable. 

We have also pointed out that, as long as seafloor 
weathering does not depend directly on surface tempera- 
ture, a climate- weathering feedback cannot operate on a 
waterworld. This is a significant result because it would 
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left with partial ocean coverage and a clement climate. 
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modeling, it would be profound, because it implies that 
waterworlds that form in the habitable zone have a path- 
way to evolve into a planet with partial ocean coverage 
that is more resistant to changes in stellar luminosity. 
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